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^ Abstract 

Gradient catastrophe and flutter instability in the motion of vortex filament within the localized 
induction approximation are analyzed. It is shown that the origin of this phenomenon is in the 
gradient catastrophe for the dispersionless Da Rios system which describes motion of filament with 
slow varying curvature and torsion. Geometrically this catastrophe manifests as a rapid oscillation of 
a filament curve in a point that resembles the flutter of airfoils. Analytically it is the elliptic umbilic 
singularity in the terminology of the catastrophe theory. It is demonstrated that its double scaling 
regularization is governed by the Painleve-I equation. 



PACS:47.32.C, 02.30.1k, 47.35.Jk 
MSC: 76B47, 58K35, 35Q05 
Keywords: Vortex Filament Dynamics, Gradient catastrophe, Flutter 

> 

•£) ; 1 Introduction 
O 



in 
o 



Motion of a thin vortex filament in an incompressible inviscid fluid in an infinite three-dimensional domain 
within the localized induction approximation (LIA) is governed by the simple equation for the induced 
velocity v 

Xt(s, t) = v(s, t) = X s A X ss = K(s, t)b (1) 
which implies the following intrinsic equations for the curvature K and torsion r 



X i K t = -2K s t - Kt s 



T t = KK S - 2tt s 



KsA (2) 
K 



Here X(s,t) is a position vector for a point on the curve representing the filament, t is time, s is the 
arclenght parameter, subscripts indicate the differentiation with respect to the indicated variables, and b 
is the binormal. 

Equations (Q} and |2]) have been derived by L. S. Da Rios in 1906 pQ and have been rediscovered sixty 
years later in (2, 3 . The detailed history of Da Rios (DR) system is presented in [4] For the localized 
induction approximation see e.g. (SJ[B]. 

In 1972 Hasimoto has observed [7] that equation (JTJ) is equivalent to the focusing nonlinear Schrodinger 
equation (NLS) 

t^t+^a« + 5lVlV = (3) 

via the transformation 

Only one year before Zakharov and Shabat [8] discovered that NLS equation (j3]) is integrable by the 
inverse scattering transform (1ST) method. Hasimoto's result has demonstrated that the whole powerful 
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machinery of the 1ST method (solitons, infinite sets of integrals of motion, symmetries, etc.) is applicable 
to the vortex filament dynamics and has lead to the explosion of interest to the vortex filament dynamics. 
Since that time various aspects of this dynamics have been analysed by different methods (see e.g. [S]- 

EH). 

Stability of vortex filament motion and formation of singularities are two important problems partially 
addressed during this period (see e.g. [13l[T5l[22l[26]). Modulation instability of the focusing NLS equation 
([3]) (see e.g. {28] Ell [30]) definitely is a key element in such an analysis. Recently it was shown that the 
dispersionless limit of the NLS equation ([3]) is quite relevant to the study of the modulation instability 
[S] - [3S]- In this limit the focusing NLS equation is the elliptic system of quasilinear PDEs. The Cauchy 
problem is ill-posed for this system and solutions exhibit the gradient catastrophe type behavior at finite 
time. Behavior of solutions of the dispersionless NLS equation near the point of gradient catastrophe and 
beyond it has been studied in [36] [37] [38] too. 

In the present paper we study the geometrical and analytical implications of such an analysis for the 
vortex filament dynamics. For this purpose it is convenient to deal with Da Rios system ^ directly. We 
first consider a dispersionless (dDR) system which describes motion of vortex filament with curvature 
and torsion slowly varying along the filament and in the motion. It is the system of quasilinear equations 
of elliptic type which in terms of the complex- valued Riemann invariant j3 = —T + iK and slow variables 
x = es, y = et, e -C 1 is of the form 

& = ~(30 + ]5)j8x. (5) 

This system well approximates DR system © until the derivatives f3 x , f3 y are not large. Hodograph 
equations for equation describe critical points of the function W which obey the Euler-Poisson- 
Darboux equation E (|, |). 

Gradient catastrophe for the dDR system happens in a point (xoj2/o) f° r given initial data. At this 
point /3 X and /3 y (or K X1 K yi r x , r y ) explode. The acceleration a = v t explodes too. So the filament 
becomes fast oscillating near the point Xo at the "time " yo. Numerical analysis (borrowed from |33| ) 
shows that such oscillations begin to expand along the filament. Intrinsically, the filament begins to 
oscillate around the rectifying plane near the point xq- So, the gradient catastrophe for the system ([5]) 
gives rise to fast oscillation which can be referred as filament flutter by analogy with airfoil flutter (see 
e.g. [39]). 

Analysis of behavior near the point xq, yo of gradient catastrophe shows that j3 — /3q +e/3*, x — XQ+ex* , 
V = Vo + e < 1, and 

W = W + e 3 (x*{fU + gV) + ^U 3 -UV 2 \ (6) 

where U + iV oc j3 and /, g are some constants. Thus the function W exhibits an elliptic umbilic 
singularity behavior in the terminology of Thorn |40j . Regularization of this singularity is one of the 
issues of this paper. 

At the point of gradient catastrophe the approximation ([5]) becomes invalid and should be substituted 
by the full DR system @. First order corrections can be obtained by the double-scaling technique together 
with appropriate modification of the function W(/3, (3) to a functional W* such that the equations of 
critical points for W are substituted by the Euler-Lagrange equation for W* . The resulting equation for 
small corrections is equivalent to the Painleve-I (P-I) equation 

% = 6ft 2 - f . (7) 

The result of the paper 36 allows us to conjecture that any generic solution of the gradient catastrophe 
behaves as j3 = /3q + ef3* where the correction (3* is described by the tritronquee solution of the P-I 
equation. 

The paper is organized as follows. Dispersionless DR system, its different forms, hodograph equation 
and other its properties are discussed in section 2. Gradient catastrophe for the dDR system is analysed 
in section 3. Geometrical implications of gradient catastrophe and flutter of filament are considered 
in section 4. In section it is shown that near to the singular point filament exhibits the elliptic umbilic 
catastrophe behavior. Regularization of this singularity via the Painleve-I equation is presented in section 
6. 
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2 Dispersionless Da Rios system 



In order to describe dispersionless (or quasiclassical) limit of DR system one, in a standard manner, 
introduces slow variables x = es, y = et with K 1 and assume that curvature and torsion are smooth 
functions of slow variables K — K(x,y) and r = r{x,y). Under these assumptions equation @ take the 
form 



K y = -2K x r - Kt x 



t,, = KK X - 2tt x + e 2 



K 



Thus in the limit e — > one has the dDR system 



K y = -2K x t - Kt x 
r„ = KK„ — 2tt x . 



(8) 



(9) 



Analogously to the dispersionless NLS equation [31]-[3I] the solutions of the dDR system © well approx- 
imate solutions of the DR system in points were K x and t x are finite. As any two-component system it is 
linearizable by hodograph transformation (x, y) <H> (K, r). The characteristic velocities for dDR system 
are complex A = — 2t + iK , Riemann invariant j3 — —t + iK and in terms of /3 the dDR system is of the 
form ©• We note that solutions of the dNLS system discussed in and those of dDR system © are 
connected by the simple relations u = K 2 , v — 2r. For the geometrical consideration the system (0) is 
more convenient. 

The dDR system is known to be integrable similar to its hyperbolic version, i.e. the 1-layer Benney 
system [41]. It has an infinite set of symmetries and integrals of motion. One of the form of the dDR 
hierarchy is given by the set of equations (see e.g. [42] ) 



Pvn= — P] . n = 1,2,3,... (10) 




where p = z + pi(x, y)z 1 + p§(x, y)z 3 + Ps(%, y)z 5 + . . . is a formal Laurent series defined by the 
equations 

p 2 = (z-P)(z~p) = (z + T) 2 +K 2 , (11) 

y n are time variables and /+ denotes the polynomial part of /. The dDR system ([9]) is the first flow 
of the hierarchy (jTU]) at n = 1. All p 2 k+i{x, y) are densities of integrals of motion for the dDR system. 
They are the dispersionless limit of the densities of integrals of motion for the full DR system ^ found 
in [HI El- 
Solutions of the dDR system can be calculated via the standard hodograph equation. It was shown 
in [431 [44] that these hodograph equations are, in fact, the equations 



Wp = 0, Wj = 0, (12) 
which define the critical points of the function of the form 

W = 1 08 +P) + |(3/3 2 + 2^/3 + if) + W(p,~p) (13) 

where the function W(f3,f3) = W((3,/3) is defined by the initial data for j3 and it is such that W obeys 
the Euler-Poisson-Darboux equation E (|, |), i.e. 

2(fi-^)W^ = Wp-Wp. (14) 

Since the function W is real valued, the second equation of (fT2]) is the complex conjugated to the first 
one. Hodograph equation is 

^ = 1 + 1(3/3+^+^ = 0. (15) 
In terms of K and r the function W is 

W = -xr + y (t 2 - ^K 2 J + W(t,K), (16) 
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the hodograph equations f| 1 5j) are given by 

W K = -yK + W K 

W T = -X + 2lJT + W T 

while the Euler-Poisson-Darboux equation is of the form 

K{WkK + W TT ) + W K 

It is the axisymmetric three-dimensional Laplace equation studied by Beltrami [45] (see also [46]). Here K 
and t play the role of the radial and axial coordinate respectively in the cylindrical system of coordinates 

For the dDR hierarchy (|T0|) the function W has the form 

W ^y yi i^L z (19) 
^ Ul J 1 2ni p{z) V ' 

where 7 denotes a small circle around infinity. Euler-Poisson-Darboux equations (|14[) is quite useful in 
the study of singular sector of the 1-layer Benney and dDR hierarchies [121 SI] 



= 
0. 



= 0. 



(17) 



(18) 



3 Gradient catastrophe 

Hodograph equations (fT2]) provide us with unique solution of the dDR system if the standard condition 

Wpp Wft 

% w w 



det < wZ wT^° (20) 



is satisfied. Due to the Euler-Poisson-Darboux equations (fT4"]) on the solutions of the dDR system with 
K the function W obeys the equation — 0. Hence 

A=|^| 2 (21) 

and, consequently, the hodograph equations (|15[) are uniquely solvable if Wpp 7^ 0. In this situation the 
derivatives (3 X and f3 x are bounded together with /3. Indeed, differentiating equation (|15[) with respect to 
x, one gets the system 

Wi P W W )\P X ) + 2\1 ) °' {22) 
So for the solutions of the dDR system one has 

* - -\w, (23) 

Thus regular sector of the dDR system ([9]) is characterized by the condition [44] 

W = 0, Wpp + 0. (24) 

In this sector j3 and the derivatives j3 x , (3 y , i.e. curvature K, torsion r, and their derivatives are bounded. 

Singular sector of dDR system is composed by solutions for which Wpp = 0. For these solutions (3 
(i.e. curvature and torsion) remain bounded while their derivatives explode. Such situation is usually 
referred as gradient catastrophe (see e.g. [47, 48 ). Generic gradient catastrophe is characterized by the 
conditions [44] 

Wp = 0, Wpp=0, Wppp^O. (25) 
At K 7^ (j3 ^ P) one also has = and the conditions (1231) . in terms of curvature and torsion arc 

W K = 0, W t ={), W kk =0, W Kt =0, Wkkk^O, W ttt ^0. (26) 
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Equations W TT — 0, WkKt — 0, Wktt — 0, are consequences of (|26j) . 

Last conditions (|26|) allow us to solve the conditions Wkk = W TT = with respect to K and r. 
Substituting these expressions into the first two conditions (|26p . one gets two equations 



fx{x,y) = Q, f 2 {x,y) = 0. 



(27) 



Thus generic gradient catastrophe for the dDR system happens in a single point (xo, yo) for given W((3, 0), 
i.e. initial data for K and r [H]. For the focusing dNLS equation this fact has been first noted in [55] . 
We would like to emphasize that in contrast to the hyperbolic case where the catastrophe locus is a curve 
in x,y here the catastrophe happens in a point. 

4 Flutter of filament 

Gradient catastrophe for the dDR system means a special behavior of a filament at the point xq at time 
yo. First one observes, using formula (20) in [16 with U = W = and V = K, that the acceleration a 
for the velocity J]]) 



v t 



e{K y b + Kb y ) = e (-2K X T - Kr x )b - K x t + T 2 n 



(28) 



explodes at Xq. Thus, at the moment to = £yo sharp bump is formed at the point so = e^o of the filament. 

Local consideration reveals another peculiarity of behavior of filament around the point Xo ■ It is well 
known, that the coordinates of a curve in a neighborhood of a point in the reference system formed by 
the tangent vector t, normal n, and binormal b and origin in the point (see e.g. formula (6.7) |49j ) are 



x\ = 

X2 = 
X3 = 



~~6 



s 3 + 



K , K s „ 



(29) 



Kt 



1 

24 



(2K s t + Kt s )s 4 



Passing to the slow variable x and slow coordinates x\ — ex±, x 2 — £X2, = ex%, one gets the same 
series (|29p for slow variables with substitution s —¥ x. At ordinary point where K, r and K x , t x are 
bounded the first terms in r.h.s. of (|29[) dominate and the curve, locally, is a twisted cubic. Its projection 
on the osculating plane (spanned by t and ft), normal plane (spanned by b and ft) and rectifying plane 
(spanned by t and b) are parabola (x2 — y x i — 0): CUS P (^Kx\ — 2t 1 x\ — 0) and cubic {xz + ^-x\ = 0) 
respectively (fig. [U see e.g. |49 a ). Characteristic feature of a curve (and a filament) in an ordinary point 





Figure 1: The figure show a curve with r > 0. At r < one should reflect the curve with respect to the 
axis t in the last figure. 



(i.e., outside the points of gradient catastrophe) is that it lies always on one side of positive direction of 
the normal and on one side of osculating plane (depending on sign of torsion r). 
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Figure 2: Possible different forms of the projection on the osculating plane of a vortex filament near to 
the catastrophe point. 



At the point xo of gradient catastrophe, the behaviors of a curve changes drastically. Indeed, when 
K x and t x become large, the second terms in xi and £3 become relevant. So, parabola in the osculating 
plane may change sign or even convert into cubic curve (fig. [5]) while in the normal plane it could be a 
plane (3,4) curve and so on. So, around the point xo of gradient catastrophe a filament oscillates from 
one side of the rectifying planes to another and back. Such oscillation is quite similar to that of airfoil 
(see e.g. [39]) and one can refer to these oscillation of a filament in the point of gradient catastrophe as 
a flutter. 

At each point of a curve there is a sphere which has contact of a third order at this point. It is called 
the osculating sphere (see e.g. [4j5] )■ It has radius 

R 2 = JI S ■ (30) 



and its center has coordinates 



^o = X + ^n-^-b (31) 



where X is the position vector of a point on the curve. 

At the point of gradient catastrophe the radius of the osculating sphere and its center blow up to 
infinity. 

Oscillation of a curve around the rectifying plane and blow up of the osculating sphere are geometrical 
features of the gradient catastrophe and flutter of a filament. At the point xo,to of gradient catastrophe 
derivatives K x , t x are so large that the terms of higher order derivatives in the DR equation (J5J become 
relevant for finite e. Solutions of the full DR system ([SJ at small e are almost indistinguishable from those 
of the dDR system at t < to = eyo (for dNLS equation see [3T]-[37]). At t > to the behavior of solutions 
of dDR and DR system are completely different. Solutions of the DR system at t > to develop a zone 
of rapid oscillations which expand around the point xo- For the equivalent dNLS equation this fact has 
been observed both analytically and numerically (see e.g. |31j-|37j). 

An example of such a behaviors is given by the figure 13 of paper [33) . In the terms of K and r 
the initial data from [33] are K(x,0) = 2e~ x and t(x, 0) = — tanh(a;). Gradient catastrophe happens 
in xo = and yo — .25. At this point the corresponding solution of DR system becomes oscillate. As 
time y growth a zone of rapid oscillations expands as y/y — yo symmetrically around the point Xq. For 
the filament dynamics this effect is seen as the creation of rapid oscillations at the point of gradient 
catastrophe and their subsequent expansion along the growing in time piece of filament (see fig. 3). Such 
a behavior is an evidence of the flutter type instability of the vortex filament motion. 
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Figure 3: Typical behavior of a flutter vortex line. 



5 Elliptic umbilic catastrophe 

In order to understand better this phenomenon one should analyse in more detail the structure of the 
singularity and the regularizing mechanism for the gradient catastrophe. For the focusing dNLS equation 
such an analysis based on the e-expansion of the integrals of motions of the NLS /Toda equation has been 
performed in 36 . Here we will follow a different approach discussed recently in [50j . 

Thus we consider a neighborhood of the gradient catastrophe point xo,yo and denote the values of 
j3 at this point by /3q. Following to the double scaling limit method (see e.g. [5T]-[52]) we will look for 
solutions of the DR system near the point of gradient catastrophe of the form 



x = xo + e x 

y = yo + e y 

/3 = A> + 



(32) 



where e«l and numbers a, a, 7 should be fixed by further consideration. For the sake of simplicity we 
restrict ourselves to the case y* = 0. We first consider the function W(/3,/3) (fl~3|) . One has 



W(x + e a x*,y , A, + e"/?*) = W° + ~(/3 + P Q )x* 



+ e 7 
1 

+ T 
1 



-x + -(3p + P ) 



2"; 



;yo 



w° 00 ) p* 



3 



+ [ -x + - {po + 30 o ) + W%)0' 



r 



37 



—1.2 



— * 3 



where W° = W 



0=Po 



0PP 

, r n _ dW 
> V¥ B — r>R 



PPpl- 



W = dW 



P=Po 



PPP 

and so on 



+ 



Hodograph equation (TT5|) and conditions (f2"5"j) imply that 



x*(/3*+3/3 ) 



(33) 



Y + T (3/?0 



w° 



(34) 



At curvature K ^ the Euler-Poisson-Darboux equations (fT4f and its differential consequences 



2f3W^ + 2((3-(3)W^ = W P p 



W- 



PP' 



imply that 



W°- = 

vv pp ' 



w° - = W°— 

yv ppp ppp 



i.e. 



yo 
4 



w°- 

pp 



0, 



w° - 

v ppp 



W°— 

vv ppp 



(35) 

(36) 
(37) 
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Taking into account equations (|34|) and (|37|) . one finally gets 

W(x + e a x*,y Q ,f3 Q + e^(3*) = W + ^e a (l3 +0 o )x* + \^ + ''(P* + T)x* + \^(aP* 3 + a^) + ■■■ (38) 
where a = ^W^^. Consequently near the point Xo,yo the hodograph equations take the form 

Wp* = ^e a+ ~<x* + e 3 '<a(3* 2 = 0. (39) 

1/2 

These equations readily imply that a — 27 and p* — x* 1 . So, near to the point xq 

^-(x-xo)- 1 /*. (40) 
Thus near to the point of gradient catastrophe the function W is of the form 

W = Wo + l^Wo + + ^W* (41) 

where 

W* = ±x*(f3*+T) + l(al3* 3 + aP* 3 ) (42) 
and (3 = — r + iK. Denoting a 1 / 3 /?* = U + iV ', one gets 

W* — x*(fU + gV) + l-U 3 ~ UV 2 (43) 

o 

where / and g are certain constants. In the Thorn's catastrophe theory a function of this form is known 
as describing elliptic umbilic singularity (e.g. [3H1IS3])- So gradient catastrophe and flutter for vortex 
filament motion enter into the general scheme of Catastrophe Theory. For the focusing NLS equation the 
appearance of elliptic umbilic singularity has been observed for the first time in [36) . 



6 Double scaling regularization and Painleve-I equation 

Next step is to regularize this singularity. Within the double scaling limit method adopted also in [36] 
one begins with the full dispersive system. Performing appropriate double scaling limit one calculate 
the required regularizing terms. An approach discussed in 50; suggests to proceed in opposite direction, 
namely, to begin with the original (say dispersionless system), to modify the corresponding function W 
adding to it the differential terms of lowest order with appropriate scaling and then to require that the 
critical point equations (IT21 for W are substituted by the Euler Lagrange equations for the modified W q , 
i.e. 



5W q 

Sf3 



0. 



^ = 0. 



(44) 



Following this approach we first observe that the contributions of j3* and (3 into the perturbation W* of 
W are separated. It is quite natural to modify the function (|33[) in such a way that this separation and 

reality property W([3,f3) = W(f3, j3) are preserved. Thus, omitting inessential second term in ([4"Tj) . and 
choosing without lost of generality 7=1, one has the following modified W 



W„ = W + > 



1 



bp:« 2 +bf3 x « 



(45) 



where S and b are appropriate constants. The corresponding Euler-Lagrange equation for (3* and j3 is 

dW q ( dW q 



d/3* \d/3*„ 



0, 



i.e. 



^x* (3* + a/3* 3 



e 5 [b/3*» x »] = 0. 



(46) 
(47) 
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To balance the first term one should choose <5 = 3. 
Thus, the modified W q is of the form 



W q = W + e 3 



\x*{p + ?) + i(a/r 3 + \aT z ) + his:, 2 + 



(48) 



and the corresponding equation for /3* is 



b{3* x * x , = a/3* 2 + ix*. (49) 
This equation is converted into the classical PainleveT equation ([7]) by simple change of variables 

x* = \£, P* = ~^ (50) 

with a\ 5 = -12b 2 . 

So, the regularized behavior at the point of the gradient catastrophe of DR system, i.e. for the 
vortex filament dynamics is governed by the PainleveT equation. The relevance of this equation for 
the NLS/Toda system near to the critical point has been observed in a different way in 36 . The 
correspondence between the solution of equation (|49]) and those given by the formula (5.21) in [36] is (at 
b = 1) u = K 2 , v = It and 

~=- l e i ' 5 (re**)- 1 '* I — ^"x* 



(re**)" 1/5 (- 



36Aq 



(3D 



where x, (3, r and ip are defined in (2.19) and (2.22) of [36]. One concludes also that e = e 2 ^ 5 . In the 
paper ([36]) it was conjectured that any generic solution of he NLS equation near to the critical point 
is described by the tritronquee solution of the PainleveT equation. Extension of this observation to our 
case allows to formulate the 

Conjecture 6.1. Generic behavior of the vortex filament near the point Xo,yo of gradient catastrophe is 
given by 

where £Iq is the tritronquee solution of the PainleveT equation (0. 

Properties of the solution ([52]) as well as other features of flutter for vortex filaments will be discussed 
elsewhere. 
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